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Using the nonperturbative method, we study the exact entanglement dynamics of two two-level dipole-dipole 
interacting atoms coupled to a common non-Markovian reservoir with different coupling strengths. Besides 
analyzing the conditions for the existence of steady state, we find, though the dipole-dipole interaction could 
destroy the stationary entanglement, that relatively strong dipole-dipole interaction does suppress the disentan- 
glement in the initial period of time. These results are helpful for the practical engineering of the entanglement 
in the future. 
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I. INTRODUCTION 

Quantum entanglement is highly relevant to the fundamen- 
tal issues of quantum mechanics and plays a central role in the 
application of quantum information science. In recent years, 
the practical schemes for the creation of quantum entangle- 
ment have gained significant improvements both in theoretical 
research and in experimental realization. One of the most pop- 
ular schemes involves a couple of spontaneously emittingtwo- 
level atoms (qubits) trapped in a single-mode cavity 1^, |3Q, 
which describes both atoms or ions trapped in an electromag- 
netic cavity and circuit-QED setups. The corresponding in- 
vestigations reach notable results on the dynamical evolution 
of the bipartite entanglement d, S and the steady states 
of the composite two-atom system 12l|30,0,[^]. 

Most of the previous studies base on an essential assump- 
tion that the coupling strengths of the two atoms to the cavity 
mode are equal [1,4], which is a good approximation for the 
atoms in microwave cavity. However, in optical cavity, the 
two atoms usually no longer possess the equivalent coupling 
strengths to the cavity mode JH 01 ■ In 0, the case of dif- 
ferent coupling constant is considered, but the dipole-dipole 
interaction is neglected in the very beginning, as in J7[] . 

In this paper, we study the two two-level dipole-dipole 
interacting atoms in an optical cavity via a nonperturbative 
method Hi, OH, 

and focus on the asymptotic and transient 
entanglement dynamics of the composite two-atom system. 

For stationary entanglement, when the dipole-dipole inter- 
action exists and the coupling strengths of the two atoms to 
the cavity mode are different, the final state of the two-atom 
system will become completely separable and the bipartite en- 
tanglement will collapse to zero, which indicates that there is 
no decoherence-free subspace in this case 111 ill . This result 
holds under both the Makovian and non-Markovian treatment. 
Thus, we conclude that the memory effect resulted from the fi- 
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nite correlation time of the reservoir in non-Markovnian limit 
does not influence the asymptotic bipartite entanglement. 

It is shown further, for transient entanglement, that strong 
dipole-dipole interaction can suppress the decoherence and 
disentanglement. However, when the strength of the dipole- 
dipole interaction and coupling strength of the atoms to the 
environment, which determines the memory effect, are in the 
same magnitude order, the dipole-dipole interaction accel- 
erates the decoherence and disentanglement for some initial 
states. 

Finally, by transforming the true mode Hamiltonian which 
is used at the start |[TolfT2ll to the quasimode Hamiltonian ifioll 
which separates the memory effect and the damping effect of 
the reservoir explicitly, we present a concrete physical inter- 
pretation for all the results derived in this paper. 

II. THEORETICAL FRAMEWORK 

Consider two two-level atoms with a common zero- 
temperature bosonic reservoir, and take their dipole-dipole in- 
teraction into account from the very beginning. Under the 
rotating wave approximation (RWA), the Hamiltonian of the 
composite two-atom system plus the reservoir is given by 

(» = i)fi3lii: 

H — Hq + H in , + Hcid, 

with 

Ho — a)icr+ (r_ + ojicr^cr^ + I dcok^kb' (u)k)b(a>k), (1) 



Hint = (aicr\> + ar 2 oV) dto k g(a> k )b(co k ) + h.c, (2) 

xj — DO 

Hdd - CT- + CT + cr_). (i) 

Now, we use the true mode Hamiltonian. Here, cr± and 
u)j are the inversion operators and transition frequency of the 
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jth qubit (j = 1,2), and b{cL>k), b'icok) are the annihilation 
and creation operators of the field mode of the reservoir. The 
mode index k contains several variables which are two or- 
thogonal polarization indexes and the propagation vector k. 
To measure the coupling strength of the atoms to the cav- 
ity mode determined by the atom's relative position in the 
cavity, we introduce the dimensionless constant aj 0]. The 
static part of the dipole-dipole interaction is proportional to 
K = [d-d-3(d- r\i)(d ■ rn)/rj 2 ]r^, where r\i - n~r 2 is the 
relative position, c/is the electric dipole moment of the atom. 
For an initial state of the form 



|-A(0)> = (cio k>i \g) z + C20 \g\ \e) 2 ) |0> £ 



(4) 



since [H,N] 

(1) (1) , (2) (2) 

cr . cr _ + cr .'crj 



0, where N = J du k W (tOk)b{LOk) + 

the time evolution of the total sys- 
tem is confined to the subspace spanned by the bases 
fl«>il*>2lO>i?,l«>il«>2lO>i,l*>il»>2ll*>£}: 

\m) = c I (fK'' W0 'k)ilg) 2 |0) £ + c 2 (fK , ' W0 '|g>ik> 2 |0> £ 

/-»CO 

+ dukC^e-^'lgWg^lWz, (5) 

\J — co 

where | ljfc>£: is the state of the reservoir with only one exci- 
ton in the A:th mode. Here, we consider the case in which the 
two atoms have the same frequency, i.e., cj\ - co 2 - wo- Ac- 
cording to the Schrodinger's equation, the equations for the 
probability amplitudes take the form 



/ci(f) = a\ 



f 



dco k . 



-i(w k -ti>o)t 



g{co k )c Wk (t) + Kc 2 (t), (6) 



ic 2 (t)=a 2 I dcuke-^-^'gito^c^ + Kdit), (7) 



where the responding function takes the form: 

co 

Suppose the atoms interacting resonantly with the reservoir 
with Lorentzian spectral density 

J(co k ) = \g(co k )\ 2 = W 2 A/n[(co k - coo) 2 + A 2 ], 
by employing Fourier transform and residue theorem, we get 
the explicit form f(t - t{) = W 2 e~ x(f ~ h \ where the quantity 
1 /A is the reservoir correlation time. 

In fact, we can obtain the solution by Laplace approach di- 
rectly. However, in order to construct a lucid physical insight, 
we choose the pseudo-mode approach J71 |9j llOfl. The corre- 
sponding equations are as follows: 



ci(f) = -iWaib(t) - iKc 2 (t), 



c 2 (t) = -iWa 2 b(t) - iKciif), 



(11) 



(12) 



b(t) = -Ab(t) - iW[aiCi(h) + a 2 c 2 (h)], (13) 



where 



bit) 



Jo 



dhe- A{ '- h) W[aic x {tx) + a 2 c 2 (h)l 



Taking Laplace transform, C\ 2 {t) '-. C\ 2 {s), b(t) = b(s), we 
get the resulting equations 



iCa> k (i) = laici(t) + a 2 c 2 (i)]g*(co k )e 



i(iOk-aM})t 



ci(t) 



c 2 (t) 



f dt K 

Jo 



f(t - t\)a\\a\c\(t\) + a 2 c 2 (h)] 



- iKc 2 (t), 



I dtif(t - fi)or 2 [o;iCi(?i) + a 2 c 2 (h)] 
Jo 

iKci(t), 



(8) 



Eliminating the coefficients c Wt (f) by integrating Eq. dHJ 
and substituting the result into Eqs. (O and (|7), we get: 



(9) 



(10) 



sci(s) - cio = — iWa\b{s) - iKc 2 (s), (14) 



sc 2 (s) - C20 = -iWa 2 b(s) - iKa(s), (15) 



sb(s) - b(0) = -Ab(s) - iW[aiCi(s) + a 2 c 2 (s)]. (16) 



Further, we introduce vacuum Rabi frequency R = W(a 2 + 
o^) 1 ^ 2 and relative coupling strengthes rj = aj(a 2 + a 2 ) ~ 1 ^ 2 
{j = 1,2). Applying the inverse Laplace transform, we get the 
solutions with the initial condition b(0) = 0: 



ci(t) = ^VaiUs-sd 



c w [s(s + A) + R 2 r 2 } - c 20 [iK(s + A) + R 2 r x r 2 \ 
s 2 (s + A)+ R 2 s + K 2 (s + A)- 2iKR 2 rir 2 



(17) 
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,. | c 20 [s(s + A) + R 2 r 2 ] - c w [iK(s + A) + R 2 ri r 2 ]) 

C2 « = L J™ <' ~ *> ,2 (5 + /t)+ ^ + ^ + /t) _ 2/ ^ 2r , r? k • (18) 



2 (s + A) + R 2 s + K 2 (s + A) - 2iKR 2 nr 2 



, -dVt I \ ( r i c io + ^20)* - lK(r\C20 + r 2 c l0 ) I 

b(t)--iR) \\m {(s - sd^, z r ^e 1 ''. (19) 

w Zj. M . 5 ,r s 2 (s + A) + fl 2 s + K 2 (s + A)- 2iKR 2 n r 2 I 



2 (s + T) + fl 2 s + K 2 (s + A)- 2iKR 2 r x r 2 



where s, belongs to the roots of the following equation for s: 



s 2 (s + A)+ R 2 s + K 2 (s + A)- 2iKR 2 nr 2 = 0. (20) 
I 



Asymptotic analysis. — In physical problems considered in 
this paper, c\ >2 (s) have poles which lie either on the imaginary 
axis, giving an oscillatory term, or in the left-hand half plane, 
giving a term which exponentially decays in magnitude. Con- 
sider function Ci >2 (s) each having a single simple pole sq = i0 
lying on the imaginary axis, and all other poles lying in the 
left-hand half plane. The residue at sq = i6 is 

lim(i - i6)c h2 (s) expOW). (21) 

s— >i6 

The other residues are of the form lim /i, exp(s,f), where hi 

is a constant determined by s,. Since each s, has a negative 
real part, the residue decays exponentially in magnitude as t 
increases. So, whether the values of c\ and c 2 are zero as 
t — > oo depends on wether Eq. (1201 has a pure imaginary 
solution. Hence, it is easy to figure out that the conditions for 
the existence of the steady state are 

K = or or, = a 2 , (22) 

which means unless there is no dipole-dipole interaction or 
the coupling strengths of the two atoms to the cavity mode are 
equal, ci(f) and c 2 (t) will tend to be zero as t — > oo. 

In the fie), \e) 2 , \e) l \g) 2 , |g>, \e) 2 , |g>, \g) 2 ) basis, the re- 
duced density matrix of the two atoms is given by: 

(0 ^ 

| Cl (f)| 2 Cl (t)c* 2 (t) 
Pa{> c 2 (t)c\{t) \c 2 {t)\ 2 

U 1 - | Cl (0| 2 - |c 2 (0| 2 ) 

(23) 

We choose the concurrence C(t) II 1 511 to measure the entan- 
glement of the two atoms. According to the definition, the 
expression of the concurrence of p a is given by lfl5ll : 

C(f) = 2|ci(0c£(f)| =2|ci(f)l|c2(f)|. (24) 

First, consider a\ = a 2 , that is, the coupling constants of the 
two atoms to the cavity mode are the same. The asymptotic 
entanglement evaluated by the concurrence is given by: 

C(t -» oo) = |c 10 - c 20 | 2 12. 

This indicates that stationary entanglement is only determined 
by the initial state of the system. 



Next, consider K — 0, that is, the dipole-dipole interaction 
is neglected. The asymptotic entanglement is given by 

C(t -> oo) = 2 \a^c w - c 20 \ \a R c 20 - c w \ /(«« + a R 2 )> 

where or = r\jr 2 . This indicates that the entanglement is 
determined by the initial state and the ratio of the coupling 
constants of the two atoms to the cavity mode. This situation 
has been discussed in details in J3]. 

Dynamical analysis. — In this part, we will focus on the en- 
tanglement dynamical evolution of the composite system in 
a good cavity, i.e., for % — > oo with % = R/A. In FIG. Q] 
we show the concurrence as a function of r = At in the good 
cavity limit. We compare the dynamics of two maximally en- 
tangled states (MES: |^±)) for four different values of relative 
strengths of dipole-dipole interaction, namely, 7C = 0, 2, 7, 20, 
with 7C = K/A which represents the ratio of the strength of the 
dipole-dipole interaction to the coupling strength of the atoms 
to the cavity field. The explicit expressions of \<p-) and \tp + ) 
are as follows: 

\<P±) = -Ul«>i \8h ± \8)i \e)i)- (25) 
V2 

We choose the coupling parameter r\ as V3/2, which is a 
typical constant in our case. Other values of r\ which exclude 
r\ - 1/ v2 show qualitatively similar behavior. Here, We 
note that when r\ = r 2 = 1/ V2, \ip-) spans an eigenspace 
of the Hamiltonian of the composite system with and without 
dipole-dipole interaction, which means that once the initial 
state is \(p~), both the population and concurrence will stay at 
their maximal value 1 as time increases irrespective of what 
the relative strength 7C is. 

For a MES \<p~), and for the case of 7C = 0, the population 
and concurrence oscillate at the very beginning, but as time in- 
creases the composite two-atom system collapses to a steady 
state, which means that the black plots for both the popula- 
tion and concurrence will converge to a straight line with a 
fixed value as r increases. However, when we take into ac- 
count the dipole-dipole interaction (7C + 0) with r\ + r 2 , both 
the population and concurrence perform damped oscillations 
and tend to be zero inevitably as time increases which means 
that the composite system has no steady state and will lose 
their bipartite entanglement completely just as our asymp- 
totic analysis describes. When 'K - 2, the entanglement col- 
lapses faster with the growing strength of the dipole-dipole 
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(a) (b) 




FIG. 1 : (color online) Time evolution of the concurrence 
and population, with the initial state being (a) \ip~), and 
(b) \ip + ), both with R = 10, A = 1 and n = VJ/2, for 
the cases of (i) no dipole-diploe interaction % = (solid 
curve), (ii) weak dipole-dipole interaction % = 2 (short- 
dashed curve), (iii) intermediate dipole-dipole interaction 
% = 7 (long-dashed curve), and (iv) strong dipole-dipole 
interaction % = 20 (long-short-dashed curve). 
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interaction. On the contrary, when "7C = 20, the entangle- 
ment collapses more slowly with the increasing strength of the 
dipole-dipole interaction. The fastest disentanglement hap- 
pens when the strength of dipole-dipole interaction and the 
coupling strength of the atoms to the cavity field are in the 
same order of magnitude (% ~ 7), given all other conditions 
the same, as shown in FIG.Q](a). Whereas, for a MES \tp+), the 
disentanglement becomes slower as the dipole-dipole interac- 
tion become stronger, which means the dipole-dipole interac- 
tion can inhibit decoherence and disentanglement, as shown in 
FIG.[T](b). In the practical engineering of entanglement, this 
result suggest that if it is hard to gain equal couplings between 
the atoms to the cavity, one can simply choose an appropriate 
distance to make the couplings of the atoms to the cavity to 
be much smaller or bigger than the dipole-dipole interactions. 
This will, at least, keep the entanglement alive for a longer 
time. 

In FIG. [1] we choose A/W = 0.1, which can be well real- 
ized within the current experimental condition llrjll . A more 
intuitive impression of the role of the dipole-dipole interaction 
can be get from FIG. [2] The explanations of these results will 
be given in details in the following section. 

III. EXPLANATIONS BY QUASIMODE HAMILTONIAN 

In order to erect a concrete explanation for these phenom- 
ena, we convert the true mode Hamiltonian with Lorentz spec- 
tral density into the quasimode form, because in the true- 
mode Hamiltonian, the coupling between the atoms and in- 
finite modes makes it hard to separate the damping effect and 
the memory effect. Using the method in II 1011 . we get: 

H = H' {) + H- m + Hdd + H 'damping, (26) 
n Q — OJ\(T + <T_ + W2<X + CT_ + OJqU a 

+ I dAAc\A)c(A), (27) 



H 'in, = W[( ai cr^ + a 2 of)a + (a l( r m + a 2 (T { ?)a\ (28) 



H damping = (A/n) l/2 I dAW c(A) + ac\A)], (29) 

\J — oo 

where c ' (A), c(A) are the creation and annihilation operators 
of the continuum quasimode of frequency A, Hdd and the other 
parameters are the same as before. 

From the above Hamiltonian, the composite two-atom sys- 
tem only interacts with one discrete mode and their coupling 
coefficients are just the transition strength a/W (i = 1,2). 
The discrete mode interacts with a set of continuum modes 
ifloll and their coupling strength only contains the width of the 
Lorentzian spectral density A which is a constant. This means 
that if we let the two atoms and the discrete mode be a new 
system and the continuum modes be the reservoir, the behav- 
ior of the new system is exactly Markovian. The physical in- 
terpretation of the quasimode Hamiltonian has been discussed 
in details in [ 10], and more details about the relations between 
the quasimode and the true mode Hamiltonian are given in the 
appendix. Now we use the quasimode Hamiltonian to explain 
the results derived above. 

The first thing is the origin of the Non-Markovian mem- 
ory effect. As we know, the memory effect originates from 
the finite lifetime of the photon in the cavity |7j, 1 1 Qfl . Accord- 
ing to Eqs. (26)-(29), the atoms only interact with the discrete 
mode, while the irreversible process only happens in the inter- 
action between the discrete mode and the continuum modes. 
So, the reabsorbing phenomenon that causes the entanglement 
to oscillate, as shown in FIG. Q] just results from the interac- 
tion between the atoms and the discrete mode, which is a re- 
versible process. Yet, due to the coupling between the discrete 
and the continuum modes, the photon escapes into the contin- 
uum modes and never comes back, for the process is exactly 
Markovian. This means that the memory effect is a finite time 
phenomenon, which fits the results in FIG.[T]where the oscil- 
lating phenomenon is only obvious at the very beginning. 
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The second thing is the disentanglement suppressing effect 
in the strong dipole-dipole interaction. From Eqs. (26)-(29), 
we can see that the atoms are only coupled to the discrete 
mode, while the damping process only happens between the 
discrete and the continuum modes. This indicates that if the 
exciton only exists in the atoms, the damping process does 
not happen. The damping of the atoms happens only when 
the exciton in the atomic system is transferred to the discrete 
mode by the interaction between them and then photon es- 
capes to the continuum modes before it comes back to the 
atomic system again. That is to say, the probability of the 
exciton transferring from the atom to the discrete mode is 
related to the disentanglement. As shown in FIG. |2] if the 
dipole-dipole interaction is much stronger than the interaction 
between the atom and the discrete quasimode, the disentan- 
glement happens very slowly. So, in terms of quantum per- 
turbation theory, the dynamics of the atoms is mainly deter- 
mined by Hdd- Because the initial state \<pJ) form an invariant 
subspace of H^, the atoms will stay in this state forever in 
the zeroth order perturbation consideration, and the disentan- 
glement is caused only by the higher-order corrections, which 
are relatively small quantities. Therefore, the disentanglement 
happens slowly. However, when the two strengths are in the 
same order of magnitude, H mt is not small quantity compared 
to Hdd, and the higher-order correction can not be treated as 
perturbation, and thus the disentanglement happens very fast. 
The analogous discussion goes well also for the case of the 
strength of dipole-dipole interaction K being much smaller 
than the strength of the interaction between the atom and the 
discrete quasimode W. 



glement, strong dipole-dipole interaction does suppress dis- 
entanglement effect under some conditions. These results are 
helpful for the practical engineering of the entanglement in 
the future. 



APPENDIX: THE RELATIONS BETWEEN THE 
QUASIMODE HAMILTONIAN AND THE TRUE MODE 
HAMILTONIAN 

The transform relation from the true mode Hamiltonian to 
the quasimode Hamiltonian is given by: 



b(oj k ) = , lA/jTy ' ., !</ + (. 



(co - ojo) + iA 

+ w*) 1/2 1. 



J — o 



to - A 



dAc'(A)]. 



(30) 



The transform relation from the quasimode Hamiltonian to the 
true mode Hamiltonian is given by: 



= Wn) 1/2 r 



dco k 



b((jL> k ) 



a> k - loq - iA ' 



(31) 



c(A) = (A/ 



doJk 



b(oj k ) 



f 



dco k 



{oj k - A)(u k -ojq- iA) 



' to k - loq - iA 



b{co), 



(32) 




FIG. 2: (color online) Time evolution of concurrence as a function 
of dipole-dipole interaction strength K, for the initial state \<p_) with 
R=W,A=l and r, = V3/2. 



IV. CONCLUSION 



Using the quasimode Hamiltonian, the initial state is given by: 
|^(0)> =ci \g) 2 \0) a ® [OJ - c 20 \eh \g)i I0>„ ® \0J. Since 

CO CO 

[N,H] = 0, where N = cr^cr ' + crfcr® + aU + 
J ^ c/Ac' "(A)c(A), we can write the state of the total system 
as 

m)) = ci(<) \eh \g) 2 \0) a ® |(U + c 2 (t) |g>, \e) 2 \0) a ® |0 W > 

CO CO 

+ &(0ls>ils>2U>„®l<U 



x: 



dcoc(co) ]g)! \g) 2 \Q) a \K). 



(33) 



According to Schrodinger's equation, we find that the equa- 
tions for c\(t), c 2 (t), and bit) are just the same as Eqs. dTTb , 
(fl2l >. and (fT3l ). This indicates that the pseudomode is just the 
discrete quasimode in the situation of Lorentzian spectral den- 
sity. 



In this paper, we extend the study of the entanglement dy- 
namics to a realistic situation where the dipole-dipole interac- 
tion is not neglected and the couplings of the atoms to the field 
are different. We find, in this situation, that there is no steady 
state entanglement. This result is valid for both the Markovian 
and non-Markovian environment. Nevertheless, though the 
dipole-dipole interaction could destroy the asymptotic entan- 
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